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Abstract 

The process pe ± — > pe^j with the photon emitted along the electron beam axis is used for luminosity 
measurement at HERA. In this paper the process is calculated including one-loop QED radiative corrections. 
In the ZEUS experiment, both the electron and the photon can be detected. Therefore both photon and 
electron spectra with and without the 7 — e coincidence are analyzed. We also calculate the process pe^ — > 
pe^ l~l + which contributes to the background in the electron tagger. 



1. Introduction 

The luminosity in the ZEUS experiment at the HERA electron-positron collider is measured using the 
inclusive bremsstrahlung process pe^ — > pe ±r y at very small photon angles (0 7 ^0.1 mrad). A photon 
spectrometer system is located on the electron beam axis at a distance of about 100 m from the interaction 
point [l[. Additional detectors are used to measure forward-going electron^ with energies between 4—9 GeV. 
The electron rate is used to calibrate the acceptance of the photon spectrometer. This study is motivated by 
augmented analysis precision which makes it necessary to go beyond the classical Bethe-Heitlcr formula [2( 
used for simulation in ZEUS so far. These simulations assume scattering on a point-like, spin-less proton. 
No proton recoil and no higher-order effects were taken into account. 

In this paper we make an analysis of the processes affecting both the photon detector and the electron 
tagger in ZEUS. Two experimental configurations arc considered, with and without coincidence of electron 
and photon. We consider the process pe ± — > pe ±r y with the photon emitted along the electron-beam 
axis. The one-loop QED radiative correction (Fig. [1} to the leptonic current is calculated. The reaction 
in which a lepton pair is produced, pe^ — > pe ± l + l~ (Fig. [2J with at least one detectable electron, is also 
considered. It does not affect the photon spectrometer but is required for the proper simulation of the 
electron tagger signals. Finally we discuss the systematic error of the luminosity measurement originating 
from the uncertainty in the cross section calculation. 

In previous theoretical estimates for these processes 0, HJ good agreement with the classical formula was 
observed. The contribution of higher-order effects appeared within the 1%-error region. However, these 
calculations were made for different experimental configurations (coincidence between electron and photon) 
and for a narrow photon energy range (8 < E 7 < 14 GeV). No generator program was published that can 
be used for detector simulations. 

This study is restricted to the QED leptonic current corrections only Weak effects are suppressed by the 
Z-boson mass in the considered momentum transfer region (Q 2 < 10 _5 GeV 2 ). The two-photon exchange 
vanishes in the forward region according to the existing estimate Q . The photons radiated by the initial or 
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+ crossed graphs 




final state proton are directed along the proton beam axis and escape the calorimeter. The other processes 
which are discussed in [(| 0] may also be neglected. 

All results in this paper are compared to the classical Bethe-Heitler formula 

dE 1 E 7 E e \E' e E e 3 J \ m^E^ 2 J w 

Here E e (E' e ) is the energy of the initial(final) electron, a is the fine structure constant and r e — a/m e is 
the classical electron radius. The approximation for the angular distribution [8[ is also checked: 

_^ e 7 / 2 ) 

d0 7 ~ ((m e /E e ) 2 + e^ 2 ' 



The paper is organized as follows: The kinematics are discussed in Sec. [2] where it is shown that kinematic 
variables may differ by up to a factor 10 30 at opposite edges of the phase space. The matrix element 
calculation is discussed in Sec. |31 Using those results a generator program has been created [§] which is 
described in Sec. |U This section also describes the precision-saving algorithm for the phase space generation 
procedure. This algorithm is necessary since the computation is complicated by large numerical cancellations 
at small scattering angles. The results are shown in Sec. [5j The paper closes with a conclusion. 

2. Kinematics 

We denote the initial and final particle momenta as follows: 

P(pi) + e ± (k 1 ) -> p( P2 ) + e ± (fc 2 ) + j(k) [+ 7 (fc')] • (3) 

The energies of the final(initial) electron and the photon are denoted by E' e (E e ) and E 1 . We consider the 
phase space region with E 1 > E° ut ~ 5 GeV and the polar angle 9 7 < 0.001. The electrons can be detected 
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+ crossed graphs for I = e 
Figure 2: The diagrams for pe^ — > pe^i+Z" process. 



if 4 < E' e < 9 GeV. The initial beam energies at HERA are Q: 

E p = 460, 575 and 920 GeV, 
E e = 27.6 GeV. 

The momentum transferred by the proton is denoted by q: 

q = P2-pi, Q 2 = -q 2 > o. 

The variation of Q 2 is very large: 



) 2 ■ 

/ rain 



„,4„,2 



3 ■ 10~ 25 GeV z , 



s 2 (E e /E^ t - 1)2 



Q 



2 

max 



A(s, to 2 , m 2 e )/s w 10 5 GeV 2 



The process of the lepton pair creation is also considered (i = e,/x): 

p(pi) + e ± (fc 1 ) ->■ ^(pa) + e ± (fc 2 ) + /+(fc 3 ) + r (fc 4 ). 

It affects the electron tagger only and is expected to be a background for the Bethe-Heitler process. 
The range of Q 2 for this process is: 
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with A = A(s, to 2 , m 2 ) = Api-k 2 — 4to 2 to 2 and <5 m = (m e + 2to;) — to 2 . 




10 10 
Q 2 , GeV 2 

Figure 3: The Q 2 -distribution for the processes; pe^ — > pe ± 7 (solid line, no cuts on E' e ), pe^ — > pe^e+e" (dotted line, 
4 < E' e <9 GeV, scaled by factor xlO 2 ), pe* -> pe ± fi+^~ (dashed-dotted line, 4 < E' e < 9 GeV, scaled by factor xlO 6 ). 

The actual distributions of the Q 2 values for the considered processes are shown in Fig. [3] The effective 
Q 2 region is much smaller than the allowed region due to the matrix element. In this paper, values of 
Q 2 > 10 ~ 3 GeV 2 are in fact neglected. 



3. Matrix element calculation 



For the one-photon exchange process the cross section is written in terms of the leptonic and hadronic 
tensors: 



spins e y spins 

The leptonic tensor is composed of the electromagnetic current vectors as: 

Lnu = E 3iJ.it = E «fc 2 e r Al [A 2 , A, ...]u fcl u fcl ef+[A 2 , A, ...] 



■\u k2 . 

The plane wave approximation is used in spite of the macroscopic interaction distance at the lowest edge of 
the Q 2 region (r - 0.1 mm at Q 2 ~ 10" 24 GeV 2 ). 

The hadronic tensor for unpolarized protons is written as: 

l^uy — — e .riflw + e i< 2 5 — , ri — w <-*mi ^2 — 4to — — . 

F 1 + t 

Here and Gm are the Sachs form factors and r = Q 2 /Am 2 ,. We use the dipole parametrization: 

GE= {i + Q 2 /K%r Gm= ^ Ge - (7) 

Here A|j = 0.71 GeV 2 and fi p is the magnetic moment of proton. The precision of the dipole formula 
is very poor for Q 2 values greater than 0.1 GeV 2 [lOj. We use it here only to check whether the proton 
size effects are visible in the considered processes. As can be seen from Fig. [3J the effective region of 
the transferred momentum is Q 2 < 10~ 5 GeV 2 for the photon- and e + e~-pair production processes. The 
dipole form factor Q differs from unity by 0.003% in that region. On the other hand, modern form factor 
formulae HI II | were mostly fitted with Q 2 > O.lGeV 2 experimental data. They have a similar (1 — Ge)- 



value at very low Q 2 . Hence we may either use the simple dipole formula ([7]) or neglect proton structure 
effects accepting the value of (1 — G E(Qmax)) ~ 0.003% as the systematic error of the calculation. The error 
of our estimate for the pe — > pe/i + /i~ process (10~ 12 < Q 2 < 1 GeV 2 ) is therefore approximately 20% due 
to the form factor uncertainty. 

The symbolic algebra program ALHEP [jjj is used for analytical computations. It creates diagrams 
using the Standard Model Lagrangian, calculates matrix elements and reduces tensor integrals over the 
virtual particles' phase space to scalar ones. Finally the matrix element is written in terms of the scalar 
products of the initial and final momentum vectors. The infrared divergences are regularized using a fixed 
photon mass parameter. The contributions of the soft radiative corrections are integrated analytically over 
the photon energy in the region mi^ < E 1 < E^ H in the laboratory frame. The basics are well described 
in Ref. [13[. The E^ H energy cutoff separates the soft and hard bremsstrahlung contributions. We vary it 
in the range of 10 — 10~ 4 GeV to control the generator consistency. The LoopTools [14] package is used 
for the numerical estimation of loop integrals. It was recompiled with real* 16 precision to avoid internal 
instability at very small Q 2 . 

The matrix elements suffer from huge numerical cancellations in the considered phase space region. The 
first problem appears when calculating momentum couplings by contracting the corresponding vectors. It 
is solved as described in Sec. SJ Another problem has to do with the cancellation of peaking terms in the 
matrix element. Let us consider typical values of the scalar products for the Born process: 

Q 2 ~ p v q - p 2 -q ~10- 22 ...10~ 5 GeV 2 , 
Ifc.qMfc^Mfe-gl-lO- 6 GeV 2 , 
k r k 2 - fc-fci - k-k 2 ~10~ 6 GeV 2 , 
Pi-kx - pvk 2 - p v k - 10 4 ...10 5 GeV 2 . 
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The major cancellation appears due to pi-kj couplings. The p 2 -kj and pi-k values are easily expressed via 
the others. But the two remaining couplings (prki and prk 2 ) are independent and may cancel numerically. 
We use the substitution: 

Si = k-ki p\-k\ — k ■ k 2 pi ■ k 2 ■ (8) 

Typical values for 8\ are: 

5t ~ lO -1 - 10 _1 « ICT 9 GeV 4 . 

Thus calculating Si precisely as 5\ — ( k-k\ k\ a — k-k 2 k 2a )-p" one can avoid one of the two p\-ki couplings. 
In this way the Born differential cross section can be written as: 



da = e 6 v4iFx(Q 2 ) - A 2 F 2 (Q 2 ) 
dT 

Ax = 
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Here the remaining fci -pi -terms appear coupled to the tiny Q 2 factor. 

The numerical behavior of this formula is stable enough to be used with fast 8-byte variables in the 
analysis. Alternatively this matrix element can be written in a more compact but still numerically stable 
way in terms of Levi-Civita tensor couplings Q. 

For the higher-order processes there is no simple substitution available, and the straightforward expansion 
in terms of one major term is not possible. At least 4 independent couplings may occur as major terms in 
various phase space regions. This problem can be partially solved using the Levi-Civita tensor couplings 
method [3j. The other possible solution is to split the phase space and to adjust the matrix element separately 
for every part. This method leads to enormous growth of the code and problems during the debugging step. 
In this analysis we solve the precision problem using high-precision floating point numbers [15j . The loop 
and the hard bremsstrahlung corrections ([3]) and the lepton pair creation process ([5]) are calculated with 
16-byte variables. They allow to use the same matrix element form in the entire phase space region. 

For cross-check purposes several matrix element representations (using various sets of independent vari- 
ables) are compared. The results are similar. However, numerical precision for different representations is 
quite different. Expressions which depend on the transferred momentum q-ki -couplings are much more effi- 
cient. Some representations require at least 32-byte variables for proper evaluation. The splitting parameter 
and renormalization constants are also varied to check the consistency of the generator. 

The higher-order matrix elements are not compact enough to be given in this paper. One can find all 
the formulae (converted to Mathematica format) together with ALHEP-scripts used for their derivation in 
Ref. M- 



4. Generator 



The typical adaptive Monte-Carlo generator is based on iterative phase space splitting according to the 
distribution of the integrated function (or its derivative). At the first (integration) step, the appropriate 
splitting grid is created and at the second (generation) step, the events are sampled. However, the narrow 
peaks of the function may lead to significant integration errors. Since only a finite number of events 
is generated in every cell, the peak may remain undetected. The integration in the forward region is 
complicated by the fact that the effective region of the phase space is negligible as compared to the total 
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range of the variables. For example, the FOAM 17] generator (with default parameters) being efficient for 
other purposes failed to integrate the Born cross section © in the whole phase space region (Q}. 

Cuts may also lead to integration errors. In our case several cuts must be applied at the integration step: 



Minimum energy of the detected photon, E™ m ~ 5 GeV 



• Minimum energy of the undetected photon, E^ H ~ 10 8 • • • 10 5 GeV (soft and hard bremsstrahlung 
separator) ; 

• Energy cut on the detected electron in the / + Z~-production process: 4 < E e < 9 GeV. The cross 
section with and without this cut is about 0.1 and 4.0 mb respectively. This cut needs to be applied 
at the integration level. Otherwise the event generation is slowed by a factor 40. 

Therefore a new adaptive Monte-Carlo integration program was created. It uses the photon energy 
as integration variable, thus avoiding E y cuts. The main idea is to carefully tune the sub-space splitting 
algorithm to minimize the integration time and the amount of memory used. The distribution of the 
derivative is examined to detect narrow peaks. A separate splitting algorithm is used if a function cut is 
detected inside the cell. Generated function values are re-used for further sub-cell integration wherever 
possible. The maximum number of steps of the phase space splitting is about 110 for 2—^3 processes and 
130 for 2 — > 4 processes. The statistical error of the integration is distributed according to the function value: 
Sli oc U. For cross-check purposes, the Born process is partially integrated with the FOAM [ItJ generator. 
The FOAM results are stable in the sub region Q 2 > 10~ 9 GeV 2 and are in good agreement with our results. 

Another problem is to compose the appropriate phase space reconstruction procedure. This procedure is 
used to map a hypercube point 6 (0, 1)^ produced by the generator to the phase space of final particles. 
The standard method of factorizing the phase space to 2-body decay sub spaces [HI is used: 

dT n (p n ) = (2iry 1 dp 2 l _ 1 dT 2 {p„ -t k n +p„_i) dT n -x(p n -i). 

However, the vectors of the 2-particle sub space are not calculated in the rest frame of the total momentum p n 
since Lorentz-transformations in the forward region lead to about 6-digit precision loss. The procedure is 
described in the following: 

We generate the recoiling proton first: 

dr » = o 9 / d( t>* I Q 4 dQ~ 2 ! <&£_i / dTn-l- (10) 
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Here <p z is the free axial rotation angle for the event, p n -i and r„_i are the momentum and the phase space 
of the remaining particles (e.g. p n -i = k + k 2 for the Born process), and the term Q 4 cancels one in the 
denominator of the matrix element. The limits for Q 2 are defined in Sec. [2] for each process. The limits on 
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For the pe — ¥ pej process the most convenient order of the remaining integration is: 

*vx=r? ^ dE : . (12) 

J J (ATT) 2 \k + k 2 \ 

B m*„ 

Here 4>k + k 2 is the angle of the plane spanned by the electron and the photon with respect to the direction 
of their sum vector k + k 2 in the lab. frame. The integration limits on the energy are defined by the general 
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formulae: 
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(13) 



For the Born process here: p = p n -i = k + fc 2 , h — fc 7 , m, = 0, (p — fc^) 2 = m 2 and the B™^ is unused. 
Finally, the limits on the photon energy are: 

E™ = maxj ^- 1 "^ ^ - bVil), ^ j . ^T" = 2~ ™° (P°-i + (14) 
For the two photon bremsstrahlung process p n -i = k + fc 2 + fc' and the phase space integral is given by: 



d(k + k 2 ) 2 f f d<j> UM dE 1 f f d<j>j; +k - 2 dE y 



7 2ir J J (4ir) 2 \k + k 2 + k>\ J 7 (4^) 2 |fc + fc 2 | 

Here the limits on E 1 are defined by Eq. (fT4)l with p n _i = fc + fc 2 and the additional condition E 1 > . 
The £J' -limits are obtained from Eq. (fT4|) by the substitutions to 2 — »• (k + fc 2 ) 2 and -E 1 ™* — > (the soft 

and hard bremsstrahlung separator). 

Similar formulae are used for the lepton pair production process (p n — i = k 2 + k^ + fc 4 ): 



dT _ /" + fc 4 ) 2 /• r dE e f f d<P k - 3+k - 4 dEt- 

J 2tt J J (^Ikz + h + hlJ J {4ir) 2 \k 3 + k 4 \ 

4m 2 E™"" E min 

L e' i — 

Here the limits on E' e are given by Eq. (|13[) with p = p n -i, nij = m e , (p — ki) 2 — (k^ + fc 4 ) 2 and the cuts on 
the final electron are the E™? n , max parameters. The integration limits for the second lepton energy Ei- are 
calculated using the same formula with p = k% + fc 4 , m, = mi, (p — ki) 2 = mf and no additional cuts for the 
-production process. For the case / = e we set E"l ax — E' e in Eq. (|13p . Alternatively, one may apply 
the detector cuts on the second electron together with the factor 1/2 to avoid double-counting of identical 
particles. 

Calculating the scalar products by contracting previously reconstructed momenta also leads to large 
numerical errors. Therefore we define products in the generator using the original revalues. 

The new code saves up to 30 decimal digits of precision when compared to the typical algorithm that uses 
momenta in the two-body rest frames [l8j . The Born process converges with the standard 8-byte floating 



point variables. Sixteen- byte floating point variables (QD package [15[) are used for the higher-order process 
calculation. 

A different set of integration variables is implemented to cross check the generator. The recoil proton 
is still reconstructed according to Eqs. (110llll[) and the sequential approach from Ref. [l8| is used for the 
remaining particles. Here the energy cuts are applied on the integrating function and the center-of-mass 
frames are used for momentum reconstruction. The result of the integration agrees with the previous one 
but the generator requires 32-byte variables (and huge computer resources) for numerical stability. 

Due to the use of high-precision floating point numbers, the performance of the generator must be 
carefully tuned. Both integration and generation steps require sufficient computer resources: up to 2Gb 
of memory and about a day of the CPU time. To optimize the generator performance we make parallel 
runs wherever possible. Thus the integration of the sub processes is performed separately and events 
are generated in parallel runs using independent sets of random numbers. We also integrate the loop- 
and factorisable soft bremsstrahlung corrections using the previously created Born-process pattern. At 
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Integration: pe~->pe~y (Born) 

a~39 . 01110 . 003mb 
Duration: < lhr 
8-byte variables 



Integration: pe~->pe~yy 

-2.90±0.04mb at E SH =1CT 6 GeV 
Duration: ~16hrs 
16-byte variables 



Integration: pe ± ->pe ± l + l 

a [e + e~] ~0 . 08mb, a [|i + |T] ~10~ 6 mb 
Duration: ~12hrs 
16-byte variables 



Integration 



Parallel runs 



Integration: pe~->pe~y + Loops 
+ Counterterms + Soft bremsstrahlung 
[Born grid is updated] 
a~36.17±0.02mb at E SH =1CT 6 GeV 
Duration: 20 parallel runs of 10 hrs 
16-byte variables 




Parallel runs 



Event Generation 

Performance: ~10 6 evts/hr 
Parallel runs are available 
16-byte variables 



Figure 4: The generator program structure. The cross section values correspond to E^ ut = 10 GeV for the photoproduction 
processes and 4 < E e < 9 GeV for the lepton pair production process. 



every point in phase space the Born matrix element is within 10 — 30% of the corrected one, |Ms orn | 2 
\MBom\ 2 + 2%t(Mg orn Mi oop ) + |M^j t | 2 . Therefore we use the Born process integration grid and simply 
update the integral value in every sub cell (using the distributed computation again). The scheme of 
optimal usage of the generator is outlined in Fig. |4l 



5. Results 

We consider the process pe^- — > pe ±r y with at least one high energy photon in the final state (E > E" Lt ). 
The value of E° ut is set to 10 GeV. We also simulate the energy acceptance of the photon spectrometer [1[ 
using an approximate formula. For the simulation of the lepton-photon coincidence, an additional cut on 
the lepton 4 < E e < 9 GeV is used. The background lepton pair creation process ([5]) is also considered. 
When the two identical particles (e or 7) appear in the detectable region only the hardest one is observed. 
The PDG'09 values are used for physical constants. The results for the cross sections are given in Tabled] 
for different beam energies. All other results are for E p — 920 GeV. 

The uncertainty of the the Born process integration is less than 0.01%. The generator error for the process 
with the one-loop radiative correction (RC) is about 0.2%. This uncertainty is dominated by the precision 
of the sub process pe — ¥ pe"f"f in the neighborhood of the soft and hard bremsstrahlung separation border 
Ey ~ E^ H . The large peak of the radiative process here results in the systematic error of ~ 0.1 — 0.2% in 
the total cross section. Better precision can be achieved by increasing the integration time and computer 
resources used. However, since the total error of the ZEUS luminosity measurement is about 1 — 2% the 
theoretical error of about 0.2% is acceptable for the cross section calculation. The uncertainty for e + e~-pair 
production is about 0.1% of the total electron tagger signal. 

The Born cross section and the photon energy spectrum coincide with the classical Bethe-Heitler for- 
mula ([T]) predictions. When the radiative corrections are applied, the total cross section is still equal to the 
classical one within the integration error. But the energy spectrum is slightly different at very high Ey values 
(Fig. [5])- This is exactly the region of the high systematic error of the radiative sub-process integration. 
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Cuts 


Process 


ct(460 GeV), mb 


cr(575 GeV), mb 


a (920 GeV), mb 




Bcthe-Heitler formula ([1} 


37.505 


37.991 


39.012 


E y > 10 GeV 


pe^- — > pe ± 7 (Born) 


37.504 


37.990 


39.011 




pe ± pe ± 7 (+ RC) 


37.58 


38.06 


39.07 


E y > 10 GeV 
4 < E e < 9 GeV 


Bcthe-Heitler formula ([1} 
pe 1 * 1 — > pe ± 7 (Born) 
pe ± pe ± 7 (+ RC) 
pe^ — > pe^e + e~ 
pe^- — > pe^^L + /i~ 


8.35 
8.35 
8.41 
0.078 
- 1.1 • 10" 6 


8.46 
8.46 
8.53 
0.079 
- 1.1 • 10" 6 


8.69 
8.69 
8.81 
0.082 
- 1.2 • 10" 6 



Table 1: Cross sections for the processes considered here. The precision of the /i + fi -cross section is significantly suppressed 
by the proton form factor uncertainty at Q 2 ~ 1 GeV 2 (see Sec. [3t. 



However, this region is suppressed by the small acceptance of the photon spectrometer Therefore there 
is no significant difference in the energy spectrum after the detector acceptance is applied (Fig. [5]). 




E Y , GeV Ey , GeV 



Figure 5: The photon energy spectrum (left) and its ratio to the Bethe-Heitler formula Q prediction (right). The dashed 
line on the left plot is the shape of the classical distribution. The ratio of the Born process cross section to the Bethe-Heitler 
formula is also presented (dotted line). The discrepancy of the total cross section is within 0.2% value. 

The angular spectrum of the photons differs from the expected classical shape ^ (see Fig. [7]). However, 
the divergence of the electron beam in ZEUS (~ 0.09 — 0.21 mrad [l|) is much greater than the mean photon 
emission angle. Thus the actual angular spectrum shape should not affect the luminosity measurement. 
To estimate the expected effect we consider the size of the photon detector. We estimate the fraction 
of events outside the spectrometer depending on the angular distribution used. The exit window of the 
ZEUS spectrometer system has a diameter of about 10 cm and is located at 92.5 m from the interaction 
point 1]. To simulate the angular spread of the incoming electron beam we used a Gaussian distribution 
with a x = 21 • 10~ 4 ,<7j, = 9 • 10~ 4 . The results, ignoring the angular acceptance of the spectrometer, are 
given in Tabic [21 The effect of the photon angular spectrum is less than 0.05%. This uncertainty is below 
the required generator error and can be neglected. If one ignores the scattering angle and only considers 
7 = for all photons the difference is about 0.4%. 

The electron energy spectrum obtained from the generator is slightly above the classical formula (see 
Fig. E]). This arises from the additional e + e _ -pair production process ©■ The difference increases in the 
low--E e region up to a few percent. This contribution may affect the calibration procedure of the ZEUS 
photon spectrometer. 

The photon energy spectrum for e — 7 coincidence is presented in Fig. [9] The deviation from the classical 
formula at the edges of the plot is due to the additional momentum carried away by the undetected photon. 
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Figure 6: The photon energy spectrum (solid line) and Bethe-Heitler formula values (dotted line) after the spectrometer 
acceptance is applied. The difference is not visible. The total cross sections agree within 0.2%. 




Figure 7: The photon angular spectrum (solid line) and the classical formula J2{ values (dashed line). 



6. Conclusions 

In this paper we have considered the processes that affect the luminosity measurement in the ZEUS ex- 
periment. The basic process pe^ — > pe ± 7 has been calculated including one-loop QED radiative corrections. 
The contribution of the lepton pair production process has also been estimated. 

The following conclusions can be drawn: 

• The Born cross section and the photon energy spectrum coincide with the classical Bethe-Heitler 
formula values ((TJ) within 0.01%; 

• The angular spectrum is significantly different from the classical estimate ([2]). However, for the ZEUS 
photon spectrometer the angular effect is less than 0.05%; 



Detector radius, mm 


7 by generator 


9 7 by Eq. © 


e 7 = 


40 


5.07 


4.97 


4.46 


50 


1.44 


1.39 


1.12 


60 


0.39 


0.36 


0.23 


70 


0.14 


0.12 


0.04 



Table 2: Fractions of events (in %) outside the photon detector acceptance depending on the angular distribution used and 
the size of the ZEUS spectrometer exit window. The angular spread of initial electron beam is simulated with a x = 21 • 10 — 4 , 
cr v = 9 ■ 10" 4 . 
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5 6 7 8 9 

E e , GeV 

Figure 8: The electron energy spectrum (left) and its ratio to the Bethe-Heitler formula Q (right). The ratio without 
the e + e _ -production process is also presented (dotted line). The dashed line on the left plot is the shape of the classical 
distribution. 
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Figure 9: The photon energy spectrum (the spectrometer acceptance is applied) if the electron is detected in the region 
4 < E e < 9 GeV. 



• The contribution of the higher-order effects to the total cross section of the inclusive pe — » pe 7 
process is within the 0.2% value; 

• The rate of the electrons in the detector is about 1 — 2% above the the Bethe-Heitler formula prediction 
due to the additional e + e~-pair production process. 

We have achieved an error of 0.2% for the cross section calculation. This value may be treated as the 
theoretical uncertainty of the ZEUS luminosity measurement procedure. Further increasing the precision is a 
difficult task that requires additional higher-order effects to be estimated. However, the error of 0.2% is much 
smaller than the systematic error of the ZEUS photon spectrometer and is acceptable for the experiment. 
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